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We analyze the qualitative features in the transition spectra of a model superconductor with multiple 
energy gaps, using a simple extension of the Mattis-Bardeen expression for probes with case I and 
case II coherence factors. At temperature T — 0, the far infrared absorption edge is, as expected, 
determined by the smallest gap. However, the large thermal background may mask this edge at 
finite temperatures and instead the secondary absorption edges found at Aj + Aj may become most 
prominent. At finite T, if certain interband matrix elements are large, there may also be absorption 
peaks at the gap difference frequencies |Aj — Aj\/h. We discuss the effect of sample quality on the 
measured spectra and the possible relation of these predictions to the recent infrared absorption 
measurement on MgE>2 . 



Among the unusual features of superconductivity in 
MgB2 is the widely discussed possibility that there may 
exist several distinct energy gaps. The larger gap or gaps 
are thought to be associated with electrons in the boron 
p xy bands, frequently referred to as "2D-bands" because 
of their very weak dependence on k z , the component of k 
parallel to the c axis. These bands appear to be strongly 
coupled to a particular [E^g) phonon mode, and are be- 
lieved primarily responsible for the high superconducting 
T c (~ 40K). Interactions between these 2D electrons and 
those in other bands (specifically, the "3D" bands formed 
by the B p z states) lead to superconducting gaps in the 
spectra of these other electrons. A multiple-gap structure 
can be described either with a simple BCS framework 
|EJ, or by more detailed ab initio calculations using the 
Eliashberg theory M . The result of such theories is that 
the average 2D gap, A2D, is about three times larger than 
A3D H. These predictions are supported experimentally 
by heat capacity tunneling |j]-|9j, photoemission 

JTo| ], and penetration depth jll]] measurements. 

In this work, we present a simple model calculation of 
the transition spectra in a superconductor with multiple 
gaps, using a natural extension of the Mattis-Bardeen 
(MB) formulae We find that this model is not only 
consistent with the observed infrared absorption edge in 
MgB2 , which occurs at anomalously low frequencies rel- 
ative to the single-gap BCS prediction, but also predicts 
that a characteristic additional structure in the absorp- 
tion may be observed under certain conditions. 

In our calculations, we consider only the effects of mul- 
tiple gaps on coherence phenomena, and on the super- 
conducting density of states, and neglect nonlocal elec- 
trodynamic effects. Such an approximation is known to 
become exact in either the dirty limit (£ <C £o> where I 
is the quasiparticle mean-free path and £0 is the coher- 
ence length), or the "extreme anomalous" limit (fc£ 3> 1, 
where k is the wave vector of the perturbation). In fact, 
neither of these limits may actually be applicable to sam- 
ples of MgB2 with multiple gaps (in particular, a single 



gap would normally be expected in the dirty limit [Ell). 
Nevertheless, previous calculations suggest that the MB 
formalism gives at least a reasonable qualitative descrip- 
tion of experiment in single-gap superconductors, even 
when used beyond its nominal validity limit. The calcu- 
lation is also very easy. Thus, it is reasonable to do such 
a calculation as a first step in understanding the absorp- 
tion spectra in MgB2 . Various further corrections should 
be included in the strong-coupling limit [ fl5|| , although 
they may not change the calculated absorption spectra 
qualitatively |ll|. 

We adopt a BCS-like multiple gap model proposed by 
Liu ei aZ . In this model, the self-consistent gap equa- 
tion is written as 

a, = £«,aJ W mh( g|f )/2) <g. (i) 

where (3 = with T the temperature, Aj is one 

of the n gaps, U%j is the corresponding n x n effec- 
tive phonon-mediated electron-electron interaction ma- 
trix calculated in jjj, iVj(£) is the normal density of states 
for the ith band, and = y/Af +£ 2 , £ = e - fx 

(where e is the single-particle energy and fi is the Fermi 
energy). Finally, u>d is a cutoff energy, which is assumed 
to be the same for all n bands. The solution of the n 
equations ([!]) gives Aj(T). In a single-gap BCS super- 
conductor u>d is of order the Debye frequency. In MgB2 
the value ujd ~ 7.5meV needed to produce T c ss 40K 
is much smaller than the physically relevant logarithmi- 
cally averaged phonon frequency ui\ a — 56.2meV |p); this 
discrepancy would probably be removed by the inclusion 
of strong-coupling corrections omitted from eq. (Q). In 
any case, the BCS model predicts gaps in MgB2 whose 
ratio and temperature dependence agree fairly well with 
detailed calculations employing the Eliashberg theory of 
superconductivity 0] (cf. Fig. 1(a) below). We therefore 
use the model of eq. (1) to calculate absorption spectra 
in MgB2 , in the hope that the results will apply at least 
qualitatively. 
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To calculate the absorption coefficients, we use the 
canonical transformation approach as described, for ex- 
ample, in Tinkham p8[ . If there are multiple gaps, one 
may generalize the standard single-gap expression for the 
coherence factors (uv 1 ±vu') 2 and (uv! ^fvv') 2 by replac- 
ing A 2 in these expressions by AA' (in the notation of 
||), so that, e. g. 



(uu' =F vv) 



J\2 



2 V EE' 
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~EE' 



(2) 



Here the upper and lower signs correspond to the so- 
called case I and case II coherence factors determined by 
the time- reversal symmetry of the matrix elements. If the 
Nj(£)'s near the Fermi level are constant for each band, 
and equal to Nj(Q), and if the electron-phonon interac- 
tion is linear (both in £ and in the ionic displacements 
u), so that the upper and lower limits of integration in 
eq. (|l]) are equal, then the terms linear in £ and £' cancel 
when the integral over the coherence factors is carried 
out. Both of these assumptions may fail to some extent 
in MgB2 (in particular, the electron-phonon interaction 
in MgB2 is exceptionally nonlinear in u GIQ). We will 
nevertheless assume that the cancellation is almost com- 
plete and neglect the terms linear in £ and £' in eq. (^) 
pOfl . We also assume that the matrix elements My for 
a one-electron transition between an electronic state in 
band i and a state in band j are the same for all states 
in given bands i and j. The transition rate in the super- 
conducting state induced by a perturbation of frequency 
u> is then found to be proportional to 



a s 



E 



M tJ \ z N t (0)N 3 (0)h 



\E(E + hcu) =f AjAj-| 



(3) 



(4) 



y/EP - A 2 ^(E + M 2 

x [f(E) - f(E + hu)] dE. 

Here f(E) = l/\e E ^ kBT +l] is the Fermi function, and the 
integration in extends from — oo to +oo, except for the 
regions where the argument of either square root becomes 
negative. The absorption coefficient in the normal state, 
a at, is given by the same formula but with all Aj = 0. 

In the single-gap case, the ratio as/otN calculated 
from eqs. (3) and (4) is usually referred to as the 
Mattis-Bardeen formula. The corresponding ratio for 
the multiple-gap case is readily calculated using eqs. 
(3) and (4). The normal state value is just ajv = 
£V. |Myf iV;(0)iV,-(0)/ia;, and hence 



1 i'j 
as 

OtN 



EtiELi 1^1^(0)^(0)^' 



(5) 



with Iij given by (EJ). 

The required integrals take very simple forms at T = 0. 
In this case, f(E) = 1 for E < and f(E) — otherwise, 
and eq. (Eh becomes 



tT=0 _ 

ij 



-A* 
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Introducing k. 
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(Tllu) 2 — (Aj — Aj) 2 , we can write this integral in terms 
of complete elliptic integrals E and K . For the case I 
coherence factors, ([j]) simply reduces to 



k^E 



whereas for the case II factors 
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[the step function 9{x) appears because the upper inte- 
gration limit in (^) should always be larger than the lower 
one]. 

In order to demonstrate the qualitative features of this 
model, we have carried out the numerical integration 
for the two-gap model of MgB2 ■ Even using this 

model, the matrix elements My still remain to be cal- 
culated. Now in MgB 2 , these two gaps are thought to 
come from two disconnected Fermi surfaces. Since these 
are disjoint, the normal-state resistivity is likely to be 
determined primarily by intraband scattering. If so, the 
off- diagonal matrix elements My would be very small. 
Various scattering processes (e. g. impurity-scattering, 
electron-phonon scattering involving a large- wave- vector 
phonon) could, however, produce non-zero off-diagonal 
matrix elements in principle. We have therefore carried 
out two model calculations, based on eq. ([| In the first, 
we have simply made the crude assumption that all the 
factors MijNi(0)Nj(0) are equal, so that the Eq. (|) be- 
comes simply 



as 



Ey hj 
n 2 huj 



(9) 



where n is the number of bands. In the second, we have 
taken My = for i ^ j, as suggested by the above argu- 
ment, but still assuming all diagonal elements MaNi(Q) 2 
to be equal, so that 



as 



rituo 



(10) 



In Figs, l(b-e), we show the transition rates as cal- 
culated from these expressions, for processes with case 
I (dashed lines) and case II (solid lines) coherence fac- 
tors, and making either of the assumptions (||) or (|Io| ) 
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at T = 0.05T C and T = 0.5T C , as indicated. The case 
II coherence factors correspond to ordinary far-infrared 
absorption, or equivalently, the real part of the frequency- 
dependent conductivity. We consider two temperatures: 
T ~ 2 and 20K. At very low temperatures, the ab- 
sorption spectra for assumption (|9|) resemble a sum of 
three single-gap spectra which become nonzero at the 
three different possible values of Aj + Aj. If the ma- 
trix elements are diagonal, the spectra would resemble 
the sum of two single-gap spectra. In either case, the 
lowest frequency where absorption occurs is determined 
by the smallest gap, which we denote Aj , and occurs at 
frequency uj — 2Aj /7i. Other absorption edges, where 
ots/otN has a slope discontinuity at T = 0, occur at other 
values of Aj + Aj, but are not very prominent, especially 
for the case II spectra. For the present model, the smaller 
gap A3C « 2.5meV, and the corresponding absorption 
edge is at 2A 3 £>//cbT c « 1.4, well below the BCS value. 
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FIG. 1. (a) The temperature dependence of the two gaps 
for the model discussed in the text; (b) and (c) The frequency 
dependence of the case I and case II absorption coefficients 
for the same model (normalized to the normal state values) 
at T = 0.05T C assuming (b) eq. (§) and (c) eq. @ for the 
transition matrix elements; (d) and (e) same as (b) and (c) 
but with T = 0.5T C . 



The case II absorption coefficients at finite tempera- 
tures [Fig. 1(c) and 1(d)] show features not present in 
the lowest temperature plots, whether the transition ma- 
trix elements are assumed to obey eq. (^|) or eq. (JlC|) . 
First, there is a weak below-gap absorption for all fre- 
quencies Tilo < 2Aj . This absorption is also present in 
the single-gap case. But in the multiple gap supercon- 
ductor, this background absorption, together with the 
peak feature discussed below, may mask the case II ab- 
sorption edge at 2 A; . (This edge is still visible for case 
I spectra, even at finite T.) Instead, for case II spec- 
tra, one of the secondary edges at A, + Aj > 2A io may 
become more prominent than the minimum absorption 
edge. In addition, if the nondiagonal matrix elements 
are substantial as in eq. (^), there is an extra peak at 
the frequency of the gap difference, hujij — |Aj — Aj|. 
At this frequency, the integrand in eq. (Q) has two mul- 
tiplicative singularities [square roots in the denominator 
of (0)] As a result, the as calculated from eqs. (3) and 
(4) is proportional to — ln|w — 0Jij\), with a coefficient 
which is proportional to the number of thermally excited 
quasiparticles. However, if the matrix elements satisfy 
the diagonality assumption (|l0|), this extra peak is ab- 
sent. The origin of this peak(s) is the same as that of the 
better-known peak at u> = in the single-gap case. This 
latter peak is responsible for the rise of the nuclear re- 
laxation rate 1 /Ti to a value exceeding the normal-state 
value as the superconductor is cooled through T c . As in 
the 1/Ti case, the actual height and width of the peak at 
I Aj — Aj I can be determined only when the k-dependence 
of the superconducting gaps, A,-(k), is included. 

We now move on to discuss the recent measurements 
of the infrared conductivity in MgB2 films C3] . In Fig. 2 
(adapted from Fig. 2 of Ref. [[uj), we show the measured 
real part of the optical conductivity, <7i(w), normalized 
to its value in the normal state at 40K, o"ijy(u;). These 
data show several unusual features. First, there is an 
apparent absorption edge at 5meV (well below both the 
weak- and strong-coupling single-gap BCS values), and 
a large background absorption below this edge. Second, 
the experimental data show a weaker and more slowly 
rising o~isI°~in ratio above the apparent absorption edge 
than a single-gap Mattis-Bardeen calculation (shown for 
2 A = 5 meV and T = 6 K as a dashed line). Finally, 
these spectra appear to have a characteristic structure 
between 3 and 4 meV, present only in the superconduct- 
ing state and reminiscent of the peak structure discussed 
above. However, this experimentally observed structure 
could lie within the experimental uncertainty, which is 
largest at low frequencies, where it is comparable to the 
scatter in the data points p2] |. Indeed, earlier data by 
Pronin et al p3[ is inconclusive regarding the presence 
of such 3-4 meV structure. More accurate measurements 
are needed to determine whether the peak structure is 
actually present in the MgB2 spectra. 

Clearly, the simple model of in Fig. 1 cannot be directly 



3 



applied to the measurements shown in Fig. 2. First, 
the the lod parameter of the BCS model ([!]) was cho- 
sen to give T c ~ 40 K, whereas the sample of Fig. 2 has 
T c ~ 30.5 K. At minimum, we should multiply all the 
energy and temperature parameters used to obtain Fig. 
1 by a factor of ~ 3/4, to account for this lower T c . Also, 
the factors MyiVj(0)iVj(0) need not be independent of i 
and j, as was assumed for Fig. 1 (b) and (c), and the 
factors MuNi(0) 2 need not be independent of i as as- 
sumed in Fig. 1(d); these factors could also depend on 
the polarization of the radiation. Finally the actual dis- 
tribution of gap values should be included. In the clean 
limit this distribution [A,*(k)] is predicted to be broad 
and to have several peaks in both "2D" and "3D" regions 
0. But impurity scattering should lead to "averaging" 
of the gap values in the experimental sample {I ~ 100 
A). The resulting gap distribution depends strongly on 
the relative magnitudes of various scattering rates, and 
in MgB2 these rates are thought to be small for scatter- 
ing between 2D and 3D states Jl^,Q, thus justifying the 
use of a "two-gap" model, as mentioned above. But for 
relatively small impurity concentrations, the gap values 
should still have some dispersion around the average A2D 
and A3D. The clean- limit results of Ref. suggest that 
two slightly different gaps, A2D,i and A2D,2, might form 
on the two 2D Fermi surfaces if the intraband scatter- 
ing is stronger than scattering between two different 2D 
bands, for relatively clean samples. 

To illustrate these effects, we have calculated as/ajv 
for case II coherence factors using eqs. (4) and (5), with 
parameters arbitrarily chosen so as to give a reasonable fit 
between the calculated results and the measured as/ctN- 
To illustrate possible effects of nontrivial gap value distri- 
bution, we used either two gaps (A2.D and A3X)), or three 
gaps (A2D,i, A.2D,2 and A3_d). In both cases, we have 
assumed that the off-diagonal matrix elements are sub- 
stantial PH . The calculated curves shown in the insets of 
Fig. 2 (solid line) have most of the features of experiment: 
an apparent absorption edge at a frequency well below 
the single-gap BCS value; a large background absorp- 
tion below this edge (|(| at finite T; and a more slowly 
rising as{ui) / 'evjv(w) with increasing frequency above the 
absorption edge than is predicted by the single-gap BCS 
model (inset of Fig. 2, dashed line). However, we empha- 
size that the extra sharp peaks in the insets are present 
only if the off-diagonal matrix elements My are substan- 
tial, which may well not be the case in MgB2. 

In conclusion, we have presented a simple calculation 
of far- infrared absorption in a model for MgB2, using 
a generalization of the Mattis-Bardeen formula to sev- 
eral gaps. The resulting absorption ratio ots{oj) / oin{u) 
shows a qualitative resemblance to the measured results. 
Specifically, the onset of absorption at a frequency well 
below the onset predicted by isotropic BCS theory, large 
background absorption below the apparent absorption 
edge, and a slower rise of as (oj)/aN (to) than in the single- 



gap case are all easily reproduced by the multiple-gap 
model. Thus, the recent optical conductivity measure- 
ments |l4|] appear to be consistent with the hypothesis of 
multiple superconducting gaps in MgB2 . If certain in- 
tcrband transition matrix elements are sufficiently large, 
the multiband model would lead to a peak at a frequency 
corresponding to the gap difference A2D — A313, which 
might be observable at finite temperatures in moderately 
dirty samples. 
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FIG. 2. (a) Symbols: experimental a\{u), scaled to the 
40K value in the normal state, as measured on a lOOnm-thick 
film of MgB 2 (T c « 30.5 K) [adapted from Ref. 13]. The 
dashed line shows the fit of the 6K data to the single-gap 
Mattis-Bardeen expression (with 2 A = 5 meV), as given in 
[13]. (b),(c) Calculated ai(w)/«rjv(w) for (b) a two-gap or (c) 
a three-gap model, assuming nonzero off-diagonal transition 
matrix elements as discussed in the text (solid lines) , and for 
a BCS single-gap model (dashed lines) using 2A = 5 meV. 
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